High Visibility Two-photon Interference with Classical Light 
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Two-photon interference with independent classical sources, in which superposition of two indis- 
tinguishable two-photon paths plays a key role, is of limited visibility of interference fringes with a 
maximum value of 50%. By using a random-phase grating to modulate the wavefront of a coherent 
light, we introduce superposition of multiple indistinguishable two-photon paths, which enhances 
the two-photon interference effect with a signature of visibility exceeding 50%. The result shows 
the importance of phase control in the control of high-order coherence of classical light. 



I. INTRODUCTION 

Interference is an essentially important topic in opti- 
cal physics, resulting in many interesting phenomena and 
important applications. The key physics lying behind in- 
terference is the superposition principle. After the birth 
of quantum physics, it was realized that superposition of 
multiple single-photon paths plays a key role in the tradi- 
tional optical interference phenomenon, which is usually 
known as Dirac's famous statement: "Each photon in- 
terferes only with itself. Interference between different 
photons never occurs" 

In 1956, Hanbury Brown and Twiss introduced the 
second-order correlation measurement, and reported a 
new type of interference effect between independent pho- 
tons, i.e., the bunching effect of thermal light 0. Soon 
after, it was realized that superposition of two indistin- 
guishable two-photon paths plays a key role in this type 
of interference effect [3]. As shown in Fig. [TJ for every 
pair of independent photons, there are two indistinguish- 
able paths for the pair of photons to trigger a coincidence 
count, and the amplitudes of these two paths are always 
of the same random phase cf> s i + 4>s2, where <p s i(i = 1, 2) 
is the random phase of the point source emitting photons 
si. Therefore, the interference term between the two am- 
plitudes will survive in the ensemble average, leading to 
the constructive or destructive two-photon interference. 
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FIG. 1. Two indistinguishable two-photon paths for a pair 
of independent photons SI and S2 to trigger a coincidence 
count. Dl and D2 are two single-photon detectors. 



case of two-photon interference with independent light 
sources, the coincidence counts consist of two parts: (i) 
The self-correlation part with the pair of photons from 
the same source which contributes a constant to the cor- 
relation function, (ii) The cross-correlation part with 
the pair of photons from different sources, which is dom- 
inated by the superposition of the two indistinguishable 
paths as depicted in Fig. [TJ resulting in the two-photon 
interference fringes. It is the existence of part (i) that 
will low down the visibility of two-photon interference 
fringes. For a quantum source such as the single-photon 
state, the self-correlation contribution could be elimi- 
nated, and therefore giving rise to a 100%-visibility two- 
photon interference However, for a classical light, the 
self-correlation always contributes, which makes the vis- 
ibility of two-photon interference fringes not exceeding 
50%. This property of two-photon interference with in- 
dependent random-phase sources was further discussed 
and demonstrated by Paul [f|, Ou @ and Klyshko 0. 
Nevertheless, the visibility of multi-photon interference 
fringes of classical lights based on third- and higher-order 
coherence could be higher than 50% |8l4ll|. 

In this paper, without the introduction of third- and 
higher-order correlation measurement, we explore an- 
other way to achieve high visibility two-photon interfer- 
ence with classical light. Instead of reducing or even re- 
moving the self-correlation contribution as in the case of 
quantum sources, we increase the cross-correlation con- 
tribution by introducing the superposition of multiple in- 
distinguishable two-photon paths (path number >2) to 
enhance the two-photon interference effect of a classical 
light. This could be practically realized by introducing 
a random-phase grating to modulate the wavefront of a 
coherent light. 



Later, Mandel gave a detailed theoretical analysis 
of two-photon interference between two independent 
sources with random phases, and predicted that the vis- 
ibility of two-photon interference fringes has a maximum 
value of 50% for classical light In general, for the 



' zhanggq@nankai.edu.cn 



II. RANDOM-PHASE GRATING 

The random-phase grating is shown schematically in 
Fig. [H[a). It is a transmission iV-slit mask with specially 
designed random-phase structure shown in the inset of 
Fig. [5{a). The transmission slit width is b and the dis- 
tance between neighboring slits is d. A random phase 
(n — 1)0 will be encoded on the light wave transmitting 
through the nth slit, where n is an integer and the phase 
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(j> changes randomly with time. Such a random-phase 
grating can be realized through a spatial light modu- 
lator (SLM) in practice, as we will demonstrate in the 
following. It is evident that there will be no stationary 
first-order interference for a coherent light transmitting 
through the random-phase grating since the phase dif- 
ference between every two slits changes randomly with 
time. However, the case will be totally different when 
one considers the two-photon interference. Here we will 
consider the case that a collimated coherent light trans- 
mitting through the random-phase grating, as shown in 
Fig- lib). 
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FIG. 2. (a) Schematic diagram for the designed random-phase 
grating with N slits. The inset shows the random phases en- 
coded on the light waves transmitting through the respective 
slits of the grating, in which the elementary phase <j> changes 
with time randomly, (b) Schematic diagram for detecting 
the two-photon interference of the light wave transmitting 
through the JV-slit random-phase grating in the Fraunhofer 
zone, where / is the focal length of the lens. 

Now, let's consider the two-photon paths in the scheme 
that contribute to the two-photon interference. For each 
pair of slits, the mth and nth slits, of the random- 
phase grating, there are always two indistinguishable 
two-photon paths as shown in Fig. [TJ (1) one photon 
transmitting through the mth slit goes to the detector 
Dl, while the other transmitting through the nth slit 
goes to the detector D2; and (2) one photon transmitting 
through the mth slit goes to the detector D2, while the 
other transmitting through the nth slit goes to the detec- 
tor Dl. Note that the same random phase (to + n — 2)<fi 
is encoded on the amplitudes of the two different but in- 
distinguishable two-photon paths. For a iV-slit random- 
phase grating as shown in Fig. [21 one may find that there 
are many such twin two-photon paths originated from dif- 
ferent pairs of slits (m, n), and the amplitudes of those 
twin paths with equal (to + n) will contain the same ran- 
dom phase (to + n — 2)cf>. These twin two-photon paths 
are indistinguishable in principle. In this way, we intro- 



duce multiple different but indistinguishable two-photon 
paths through the TV-slit random-phase grating, and the 
superposition of these multiple two-photon amplitudes 
would enhance the two-photon interference, leading to 
high- visibility two-photon interference for classical light. 



III. THEORETICAL RESULTS 

To clearly illustrate the two-photon interference effect 
of the light field transmitting through the random-phase 
grating, we will calculate the second-order correlation 
function of the transmitting field in the Fraunhofer zone, 
i.e., in the focal plane of a lens put behind the random- 
phase grating, as shown in Fig. (Hb). For simplicity, we 
assume that the coherent light incident on the random- 
phase grating is a plane wave and a single-mode one, in 
one dimensional case, the field operator on the detection 
plane is expressed as [H, EH 

„| jv-i 
£ (+) (^) oc / 2 e- if3 J x »dx s e" m(ftd " 0) a. (1) 

J 2 n=0 

Here a is the annihilation operator, j3j = fcsin^' and 
tan#j = Xj/f (j = 1, 2) with /, k and 6j being the focal 
length of the lens, the wave vector and the diffraction 
angle of the light wave, respectively. 

First-order Spatial Correlation Function — The first- 
order spatial correlation function in the detection plane 
can be expressed as [l2| 

G<»{x j ,x J ) = {E'(x j )E(x j )), (2) 

where E*(xj) is the eigenvalue of the field operator 
E i x j) 011 the state of source (coherent state), and (• • • ) 
represents the ensemble average. By substituting Eq. JT]) 
into Eq. ©, one gets 

Here the term sin (/3jb/2)/(/3jb/2) 2 represents the 
diffraction from a single slit. It is evident that the in- 
tensity distribution in the focal plane is a sum of the 
diffraction intensities from N different slits. There is no 
stationary first-order interference among these slits of the 
random-phase grating. 

Second- order Spatial Correlation Function - - The 
second-order spatial correlation function at the detection 
plane can be expressed as [12[ 

G^(x u x 2 ) = (E*(x 1 )E*(x 2 )E(x 1 )E(x 2 )} . (4) 

By substituting Eq. (fTJ) into Eq. (j4]), the second-order 
spatial correlation function can be deduced as 
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Note that the terms e -i^xd e -imp 2 d and e -in^d e -im^d 
correspond to the twin two-photon paths transmitting 
through the mth and nth slits of the random-phase grat- 
ing, respectively. Here we introduce the delta function 
S(m — n) to show that there is only one path when the 
two photons transmit through the same slit to trigger 
a coincidence count. Following the superposition prin- 
ciple, in Eq. ([3]), the amplitudes of all different but in- 
distinguishable two-photon paths with the same random 
phase l(p (I = 0, 1, • • ■ , 2iV — 2) are superposed to cal- 
culate their contributions to the coincidence probability, 
and then the coincidence probability contributions from 
those with different random phases l<p are added to get 
the total coincidence probability. 

Thus, the normalized second-order spatial correlation 
function can be calculated as 



c/ (2) (xi,x 2 ) 
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sin 2 ((Z' + l)(/3i -P 2 )d/2) 



(6) 



l' = -N 



sin 2 ((A - &)d/2) 



One notes that the terms sin 2 ((7' + — 
fi 2 )d/2)/ sin 2 ((/3i — /3 2 )d/2) are of the similar formula as 
the multiple-slit single-photon interference function [ljj], 
and therefore can be called as multiple-slit two-photon in- 
terference function. It is seen that g^ 2 \x\, x 2 ) in Eq. ([6]) 
is a sum of (27V— 1) multiple-slit two-photon interference 
functions sin 2 ((7' + l)(/3i - P 2 )d/2)/ sin 2 ((/?i - /3 2 )d/2) 
introduced by the random-phase grating, each one is 
associated with a group of different but indistinguishable 
two-photon paths which are characterized by the same 
random phase l(f> (I = 0, 1, • • • , 2N - 2) in Eq. ©. 
These multiple-slit two-photon interference functions 
are peaked at the same condition (f3i — f3 2 )d = ±2mr 
(n = 0, 1, 2, •••) due to the constructive interference 
effect, i.e., when the phase difference among different 
but indistinguishable two-photon paths are an integer 
multiple of 2ir. The constructive interference peak for 
each multiple-slit two-photon interference function is 
(I' + I) 2 , and therefore, one can get the interference 
peak of g^(xi,x 2 ) to be (2iV 2 + I)/(3JV), according 
to Eq. ^j. It can also be easily seen from Eq. © 
that g( 2 > (x\, x 2 ) is a periodical function of the position 
difference {x\ — x 2 ) with a period of A = Xf/d in 
the paraxial approximation, which is exactly the same 
as that of the multiple-slit single-photon interference 



pattern of a normal grating with respect to the po- 
sition x on the detection plane [lj]. On the other 
hand, the minimum of g^(xi, x 2 ) is achieved at the 
condition (ft - (3 2 )d = ±(2n + 1)tt {n = 0, 1, 2, 
• • • ) due to the destructive interference effect among 
multiple two-photon paths. However, the minimum of 
g( 2 \xi 7 x 2 ) is not zero but calculated to be 1/N due 
to the existence of the cases when the two photons 
transmit through the same slit of the grating. Therefore, 
the visibility of the two-photon interference, defined as 

V = (g£l - ff<l)/G$L + 9 { ±), is (N 2 ~1)/(N 2 + 2), 
which grows quickly with the increase of slit number iV 
and exceeds 50% when N > 2, as shown in Fig. [3J 




FIG. 3. The dependence of the two-photon interference visi- 
bility V on the slit number N of the random-phase grating. 



IV. EXPERIMENTAL DEMONSTRATION AND 
DISCUSSIONS 

Experimental Setup — Figure [5] shows the experimen- 
tal setup that we used to measure the two-photon in- 
terference effect of the light field scattering from the 
random-phase grating. In our experiments, a single mode 
laser with a wavelength of 780 nm was introduced as the 
light source, which was expanded and collimated through 
a beam expander to obtain a plane wave. The expanded 
and collimated light beam was then reflected by a beam 
splitter BS and incident normally onto a random-phase 
grating. Here the random-phase grating was composed 
of a TV-slit amplitude mask (b = 72 /jm and d = 400 /im) 
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FIG. 4. Schematic diagram of the experimental setup. A 
single mode 780-nm laser was introduced as the light source. 
A TV-slit amplitude mask (TV = 2, 3, 4 and 5, b = 72 /tm 
and d = 400 /im) and a SLM were used to construct the 
random-phase grating shown in Fig. [2] A CCD camera was 
located at the focal plane of the lens L to measure the intensity 
distribution and the second-order spatial correlation function. 



and a reflection-type phase-only SLM (HEO 1080P from 
HOLOEYE Photonics AG, Germany) which was put as 
close as possible to the mask. The light first transmit- 
ted through the TV-slit amplitude mask, and then was 
reflected back from the SLM and finally re-transmitted 
through the TV-slit amplitude mask again. The SLM pro- 
vided the desired random-phase structure on every indi- 
vidual slit of the TV-slit amplitude mask as shown in the 
inset of Fig.[Ha). At last, the light waves scattered from 
the random-phase grating were collected by a lens L with 
a focal length / = 80 cm. Both the intensity and the 
second-order spatial correlation measurements were per- 
formed on the focal plane of the lens L by using a charge 
coupled device (CCD) camera. 

Figure [5] shows the measured single-photon and two- 
photon interference patterns on the detection plane (i.e., 
the focal plane of the lens L) at different conditions, in 
which the empty circles are the experimental results while 
the red curves are the theoretical fits, respectively. 

Results for Traditional Grating — When there is no 
electric signal loaded on the SLM, our experimental con- 
figuration is essentially the same as a typical setup to 
measure the single-photon interference of a traditional 
TV-slit grating. In the experiment, we measured the 
single-photon interference patterns of the TV-slit grat- 
ings (TV = 2, 3, 4 and 5, respectively). The results 
are shown in the first column of Fig. [5l As expected, 
stationary single-photon interference fringes described 
by the multiple-slit single-photon interference function 
sin 2 (TV/3ei/2)/sin 2 (/3d/2) were observed. The period 
between the neighboring principal intensity peaks was 

I 



measured to be 1.57 mm on the detection plane, and 
(N — 2) sub-peaks appear between the two neighboring 
principal peaks of the stationary single-photon interfer- 
ence fringes. Note that the normalized second-order spa- 
tial correlation function g^ 2 \x\, x^) in this case was con- 
firmed to be a unity (not shown in Fig. [5]). 

Results for Random-phase Grating — When the SLM 
was loaded with the random phases, the random-phase 
grating was constructed. In this case, there should be 
no stationary single-photon interference fringes since the 
phase difference between every two slits changes ran- 
domly. The second column of Fig. [5] shows the exper- 
imental results. One can find that the single-photon in- 
terference fringes were erased, leaving an intensity distri- 
bution enveloped by the single-slit diffraction profile as 
described by Eq. ©. One notes that there are still some 
residual intensity fluctuations which deviate from the in- 
tensity envelop described by Eq. ([3]). This is mainly due 
to the unavoidable phase flicker effect of the SLM, which 
leads to imprecise control on the random phases on the 
slits of the grating. 

Interestingly, although the single-photon interference 
fringes disappear with the random-phase grating, two- 
photon interference fringes appear. As shown in the 
third column of Fig. El the second-order spatial corre- 
lation function g( 2 ' (xi, x 2 ) exhibits itself in the form of 
high quality interference fringes, which is in good agree- 
ment with the theoretical prediction by Eq. ([6]). Here the 
second-order correlation function is calculated through a 
formula g {2 \x u x 2 ) = (I(x 1 )I(x 2 )) / ((I( Xl )){I(x 2 ))) by 
using 10000 frames of the instantaneous intensity distri- 
butions measured by the CCD camera [T3, E3] ■ It is seen 
from the experimental results shown in the third column 
of Fig. [5] that, the two-photon interference fringes are 
peaked at the position differences x\ — x 2 = ±2mr f / (kd) 
but minimized at the position differences x\ — x 2 = 
±(2n+ l)irf/(kd) (n = 0, 1, 2, •••), respectively. The 
period of the two-photon interference fringes A was mea- 
sured to be 1.57 mm, in good agreement with the predic- 
tion of Eq. ©. On the other hand, sub-peaks typical for 
the single-photon interference fringes shown in the first 
column of Fig. [5] were not observed in the two-photon 
interference fringes in the third column of Fig. [5] This 
is due to the fact that g^ (x\, x 2 ) is a sum of (2TV — 1) 
different multiple-slit two-photon interference functions 
(see Eq. ([6])), and these different multiple-slit two-photon 
interference functions are always in phase at their princi- 
pal peaks but out of phase at the sub-peaks. Moreover, 
the visibility of the two-photon interference fringes was 
measured to be 44.9%, 59.1%, 62.3% and 71.9% for the 
TV-slit random-phase gratings with TV = 2, 3, 4 and 5, 
respectively. As predicted by Eq. ([5]), the visibility of 
the two-photon interference fringes increases with the in- 
crease of the slit number TV of the random-phase gratings 
and surpasses 50% when TV > 2. 



Single-photon Interference 



No Single-photon Interference 



Two-photon Interference 




-1 o 1 
x (mm 



-2-1 0,12 

x r x 2 (mm) 




=3 0.8 
CO 



CO 

c 

CD 



: 0.4 
0.2 
1.0 



0.8 

-5. 

f ° { 

c 


c 0.4 



-2-10123 

x (mm) 



0.2 




o Experimental Data 

Theoretical Fit g\ 


M 


2.0 - 






1.6 ■ 

CM 

X 






1.2 ■ 






0.8 ■ 
0.4 - 


-3 -2-10 1 

x (mm) 


2 3 




o Experimental Data 
Theoretical Fit n m 


M 


2.4 ■ 
2.0 ■ 






*- 1.6 - 
S" 1.2 - 






0.8 ■ 






0.4 - 


-3 -2-10 1 

x (mm) 


2 3 






3.2 
2.8 
"5,2.4 

" 1.6 
35 1.2 

0.8 
0.4 



-3-2-10123 

x r x 2 (mm) 


■ o Experimental Data 
Theoretical Fit 









x (mm) 



x r x 2 (mm) 



FIG. 5. The single-photon interference and the two-photon interference of light on the detection plane at different conditions. 
The empty circles are the experimental data, while the red solid curves are the theoretical fits. The first column shows the 
stationary single-photon interference fringes with the normal JV-slit gratings. The second-column shows the averaged intensity 
distributions over 10000 frames of instantaneous single-photon interference patterns with the iV-slit random-phase gratings. 
And the third column represents the two-photon interference fringes with the TV-slit random-phase gratings. Here we used a 
CCD camera to measure the instantaneous intensity distributions on the detection plane, as shown in Fig. [4] In all cases, the 
slit number TV of the grating was set to be 2, 3, 4 and 5, respectively. 



Further Discussions — One may note that the ran- 
dom phases encoded on the slits of the random-phase 
grating are not fully independent but indeed correlated 
with respect to each other. Therefore, our case is differ- 
ent from the case discussed by Mandel Q , where classical 
lights with fully independent random phases are consid- 
ered. For classical lights with fully independent random 
phases, the visibility of two-photon interference fringes 
cannot exceed 50%. While by appropriately controlling 



the random phase structure encoded on a coherence light 
field, one could realize two-photon interference with its 
visibility exceeding 50%. It is known that controlling 
optical phase plays a key role in the single-photon inter- 
ference effect, our results show that it may also play an 
important role in controlling the high-order coherence of 
light. 
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V. SUMMARY 



In summary, we have designed a kind of two-photon 
grating with a special random-phase structure, through 
which the single-photon interference is smeared out but 
the two-photon interference appears. With such a 
random-phase grating, superposition of multiple indis- 
tinguishable two-photon paths is introduced, which leads 
to a high visibility two-photon interference fringes of 
classical light. Theoretically, the visibility of the two- 
photon interference fringes for a coherent light trans- 
mitting through a TV-slit random-phase grating reaches 
(N 2 - 1)/(N 2 + 2). Experimentally, the visibility of the 



two-photon interference fringes with a iV-slit random- 
phase grating (N = 2, 3, 4 and 5) was measured to be 
44.9%, 59.1%, 62.3% and 71.9%, respectively. The re- 
sults show the possibility to control the high-order co- 
herence of light through optical phase. 
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